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Q ! Abstract 
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for relativistic quantum field theories using the nonperturbative Hamiltonian approach 
of Discretized Light-Cone Quantization. The method is intended for eventual applica- 
tion to quantum chromodynamics in 3+1 dimensions. Here the fundamental principles 
are illustrated concretely by treatment of QED in two dimensions. The code is in- 
Oh, tended as a basis for extensions to include more complicated gauge symmetry groups, 

! such as SU(3) color, and other quantum numbers. The code was written in Fortran 

77 and implemented on a DEC 5000-260 workstation with a typical runtime of 0.2s at 
total momentum 10. 



Program Summary 

Title of program: ms_main.f 
Catalogue number: ???? 

Program obtainable from: MPI fur Kernphysik, Heidelberg, FRG 
Licensing provisions: none 

Computer: workstation DEC 5000-260; Installation: MPI fur Kernphysik, Heidelberg, 
FRG 

Operating System: ULTRIX version 4.2 
Programming language used: Fortran 77 
Program storage required: 255kbyte 



'Present Address: Institut fur Physik/COM, Humboldt-Universitat, Invalidenstrasse 110, D-10115 
Berlin, Germany 
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Data storage required: 80kbyte 
Peripherals used: terminal 
No. of bits per word: 32 
Number of lines in program: 1894 

Keywords: particle physics, field theory, Hamiltonian method, QED, Schwinger model, 
Discretized Light-Cone Quantization 

Nature of the physical problem 

Gauge field theories are notoriously difficult to solve non-perturbatively, in particular non- 
Abelian gauge theory which underlies quantum chromo dynamics. We describe the rela- 
tively new approach of [1] to this problem, Discretized Light-Cone Quantization (DLCQ), 
which in many respects is complementary to Lattice Gauge Theory. We detail the method 
in some generality but present the programming application to Quantum Electrodynam- 
ics in 1+1 dimensions [2]. Mathematically, in DLCQ one formulates the theory in light- 
cone coordinates on a finite space-interval or volume with periodic (and/or antiperiodic) 
boundary conditions on the quantum fields. Correspondingly the plane wave Fock basis 
is discretized. A particle number and momentum cut-off lead to a completely specified 
numerical problem. One can thus solve the Schrodinger equation for bound states. The 
continuum limit must be obtained by numerical extrapolation of, for example, the invari- 
ant masses. The generalization to other theories is evident. The program is intended as 
a template for further calculations in theories approaching the complexity of QCD(3+1), 
and has been already used as such. 

Method of solution 

The algorithm can be briefly described as follows: 

Step 1: Choose a harmonic resolution and particle number cut-off. 
Step 2: Construct the Hilbert space (Fock states) respecting statistics and symmetries. 
Step 3: Evaluate for each matrix element the analytic expression, making sure all conser- 
vation laws are obeyed. 
Step 4: Add all allowed terms to obtain the Hamiltonian matrix. 
Step 5: Diagonalize, yielding the complete mass spectrum and wave functions. 
Step 6: Repeat steps 1-5 for higher resolutions to check convergence. 

Restrictions: The resolution is restricted by the available memory, as the Hamiltonian 
matrix is stored. The run time is negligibly restrictive for matrices of dimension 1500, but 
grows with matrix dimension d 3 . 

Typical run time: 0.2s for the test case of resolution 10. 
References: 

[1] H.C. Pauli, S.J. Brodsky, Phys. Rev. D32 (1985) 1993, 2001 
[2] J. Schwinger, Phys. Rev. 128 (1962) 2425 
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Long Write-Up 



1 Introduction 

The major physical theory calling for the development of nonperturbative techniques for 
solving quantum field theories is quantum chromodynamics (QCD). It is arguably the 
theory of strong interactions describing hadrons. However to extract the physics of mesons 
and baryons within a picture of relativistic bound states of quarks and gluons is extremely 
difficult. Lattice Gauge Theory Q is one nonperturbative method which is presently 
computing ratios of hadron masses with ever improving success [g] . However it offers little 
insight into how the above constituent or partonic picture emerges out of QCD. Discretized 
light-cone quantization (DLCQ) exploits the computational elegance and simplicity of the 
lattice method - by working with a momentum space lattice. However it also employs 
the physically intuitive power of Dirac's 'front form' || in which the vacuum is essentially 
trivial and one of the momentum operators is independent of the interaction and positive 
definite. Thus the low energy part of the QCD hadron spectrum in such a framework 
would be well described by a low number of light-cone Fock space quanta built from the 
simple light-cone vacuum. DLCQ obtains both the spectrum and the wavefunctions of 
relativistic bound states of the theory by direct diagonalization on the computer of a 
finite dimensional Hamiltonian matrix. How a Lagrangian field theory is reduced to such 
a problem is the whole art of DLCQ. 

This is not the place to detail the history or subtleties of DLCQ for which the reader is 
referred to Q. It suffices to say that QCD in (3+1) dimensions has not yet been attacked 
with the full power of the method. Rather it is being worked towards via progressively 
more complicated models by ourselves [|| and others || 0, This is therefore an 

opportune point in the development of the method to present the programming details 
at their simplest but with enough richness to display the power of the method. We turn 
to quantum electrodynamics in one space and one time dimension, QED(1+1), originally 
treated by Eller 0. This model shares with QCD(3+1) the feat ure of being a relativistic 



theory of gauge vector bosons coupled to fermions. It is, however, unencumbered by 
the technicalities of color and transverse dimensions. These demand a richer Fock space 
and additional Hamiltonian structures the programming of which can nevertheless be 
grafted onto the basic procedures we outline. Moreover the simplicity of the theory here 
means Hamiltonian matrix elements can be analytically calculated by hand. This may 
be technically inefficient in more elaborate theories in which case part of this task could 
also be programmed. Two aspects, however, could change the underlying method of 
programming: nontrivial renormalization and dynamical zero modes of gauge boson fields. 
These problems certainly become acute in higher dimensions and remain at the frontier of 



present research into the method [11, 12] but we will not say anything more about them 
in the present work. 
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In the following section we outline briefly how one builds from the Lagrangian defining 
the field theory to the matrix problem via DLCQ. In Section 3 there is a description of 
the various components of the program and how they can be implemented in a test run. 
The final section gives a statement of conclusions. 

2 Light-Cone Hamiltonian Matrix 

This section is intended to explain notation and concepts in the formulae and the corre- 
sponding computer codes. To show the generality of the method we describe the proce- 
dure for QCD(3+1), and later show how this concretely is worked out for the model case 
QED(1+1). 

Derivation. For a given quantum field theory in an arbitrary number of space-dimensions, 
DLCQ sets out to solve the eigenvalue equation 

H LC \V) =m 2 \^). (1) 

Here i?LC = P^Pu = 2P~ P + — P"j_ is the so-called 'light-cone Hamiltonian'. In two 
dimensions P± = 0. Evidently Eq.(|l|) is a realization of the relativistic energy-momentum 
relationship on physical eigenstates of the given theory. The (four) momentum operator P^ 
is built from the energy-momentum tensor, a conserved quantity under evolution in light- 
cone time x + = -^(x° + x 3 ). Our conventions follow Kogut and Soper jl3|| . Initial data for 
the independent Lagrangian fields are specified on the null-plane surface x + = 0. This will 
be done by expanding the independent fields at x + = in momentum modes. A momentum 
space 'lattice' is obtained by restricting the one dimensional light-cone space variable 
x~ = -^( x ° — x 3 ) to a finite interval from — L to +L and invoking boundary conditions on 
the fields. A similar demand can be made on x± in higher dimensions giving discretized 
k±. Periodic boundary conditions for bosons is the only choice that is consistent with the 
standard Euler-Lagrange equations. As fermions appear bilinearly in the Lagrangian of 
any of our candidate gauge theories, either periodic or antiperiodic boundary conditions 
can be used. The program we present in the sequel allows the user to choose. Thus the 
corresponding momenta k + are discretized in, respectively, integer or half-integer units of 
ir/L. Due to light-cone kinematics, these integers are either positive or zero. Ignoring zero 
modes of gauge bosons then, we can choose the light-cone gauge A + = 0. The Maxwell 
and Dirac equations give constraints allowing one to eliminate the gauge potential A~ 
and the lower spinor component ^l- Thus, one is brought to the independent degrees of 
freedom, which in physical dimensions are the transverse gluons A 1 , i = 1,2 and the upper 
spinor component ^r. These are subject to quantization. Quantization is achieved by 
imposing equal x + (anti)commutation relations on the dynamical fields. These translate 
into (anti)commutators on the Fourier coefficients of the aforementioned momentum mode 
expansions. The coefficients become Fock operators of boson, fermion and antifermion 
states. The above eigenvalue problem is then expressed in terms of these operators: the 
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light-cone Hamiltonian is written in second quantized form, and the wavefunctions are 
built from the Fock operators as an expansion in the number of fundamental partons with 
amplitude coefficients which are obtained by solving Eq.(|l|). The eigenvalue problem for 
the bound states is thus completely specified numerically. 

Now we specify using the model theory QED(1+1). First, in 2 dimensions there is no 
spin and, neglecting zero modes, there are no dynamical photons. Thus is the only 
independent quantum field leading to fermion and antifermion Fock operators in terms of 
which the eigenvalue problem must be specified. In the light-cone Hamiltonian one can 
separate out the dimensionful Lagrangian parameters of fermion mass rrif and coupling 
constant g by invoking the dimensionless operator functions in K, T, C, F and S 



H LC = 2K(mj + (1 - - | ; mM9 )T + - ■ ; mM9 (C + F + S) . (2) 
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1 1 



l + 7r(^)»' l + vr(^)2 



These operators are built from combinations of the fermion and antifermion Fock operators 
and analytic expressions for their matrix elements will be given in the sequel. Here K 
represents the dimensionless momentum called the harmonic resolution and T is the kinetic 
energy. In the literature, the potential energy terms C, F and S are called respectively 
'contractions', 'forks' and 'seagulls' because of a diagrammatic representation of their 
operator structure |14|, The simplicity of QED(1+1) allows application of the Wick 
expansion leading to analytic expressions for the matrix elements of i^LC between arbitrary 
Fock states. These expressions are given below. The matrix is normalised by a factor 
(mj + 

Notation. In the absence of spin and color, the only quantum number needed to specify 
states is longitudinal momentum. In the following, we use k to represent this momentum 
for single electron states and p for positron states. The quantity N will represent the 
particle number. As we will construct multiparticle Fock states, we attach a label to the 
momenta, fcj, to indicate the position of the corresponding Fock operator in what will be 
momentum ordered states. Thus, for example, k\ will always represent the electron with 
lowest momentum in the state. The summation indices of the sums thus run over 

the positions 1,2 ... N. Incoming momenta ki,Pi carry no prime while outgoing momenta, 
k'^p'i, are primed. As usual, 'incoming' means those states on the right-hand side of 
expressions. We denote incoming and outgoing states by \i) and (/| respectively. For 
convenience we surpress the forbidden zero momentum transfer i terms in the formulae. 

Because of contractions between states of like-momentum we can use the 'specta- 
tor/participant' picture for the interaction. Particles not involved in the interaction of 
up to four participants have to fulfill a pairwise momentum conservation. This will be 
represented by the symbol A s 



A «=n <5 wk{ii <5 PiP5> ( 3 ) 

* 3 

with indices i,j running only over spectator positions in the Fock space vectors. 
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The relative sign S r of graphs can be calculated using reference graphs and the 
number of transpositions of particles necessary to get to these. This leads to terms 
S r = and Cab = | x if b>a . o e i se . | . 

Matrix elements. We thus obtain the following c-number expressions for the matrix ele- 
ments of the aforementioned operators. 

(f\K\i) =5 NN > Yl 6 k n k>„ II W„ ( k n+Pn), (4) 

n=l. ..iV n=l...N n=l...N 



(f\T\i) =s NN , n s knK n ^ e Uv + ^ (5) 

n=l...N n=l...N n=l...N Z ^ n 



(6) 



</|C|i> = 6 NN , Y[ 5 knk , n Y[ 5 PnP , n x 

n=l...N n=l...N 
V 1 / 1 1 1 1 \ 

nk..N q keV 2 H*» " ^ + " (*» + (ft. + ?) 2 ' ' 

n,m,s,t=l...N 

(^ +Pm -^- P J,o(-l) ( " +m+S+ ^^^ " (^Tfc^ 

\ ^ X , ( -\\(n+m+s+t+c n , m +c s ,t) \ 

+ 2 d kn+k m -k' s -k' t ,0{-i-) _ fc /)2 

+ + D (-l) (n+m+S+ * +Cn ' m+Cs,t) 7 ^— TTo), (7) 

</|F|i) =A S ^ 

n = 1 ... AT; 
m,s,t = l...N' 

(^-^- fci - Pt ,o(-l) (w ^ + ^ + ^(- ( ^_ 1 ^ )2 ) 
+ ^-^-rt-^,o(-l) (B+m+a+f+Cm -' +Af) (- ft/ ' , )2 )) 

V/^ra Km/ 

,Ar'+2 x 

n, m, s = 1 . . . TV; 
t = 1 ... N' 

($k +k + v , fc , (-l)("+ m + s+t + c ".- +JV, )(_ 1 ) 

+^p m+ k^ P 'M- 1 ) {n+m+S+t+C ^^ (8) 

These analytic expressions were implemented in the computer code we discuss in the 
following section. 

3 Features of the Program 
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3.1 General Run Features 

The program 'ms_main.f ' was written in Fortran77 and was compiled and run on a Unix- 
system DEC 5000-260 work station, using routines of the Numerical Algorithm Group 
library NAGLIB Mark 15. It observes the standard Fortran 77 conventions, but uses the 
comment sign '!' and longer variable names which can be easily shortened if necessary 

The 'physical' input parameters are: choice of fermion boundary conditions as a logical 
constant, the harmonic resolution as a natural number, the fermion mass and the coupling 
constant as positive real numbers. The desired particle number truncation is given via two 
switches, fixed particle number and maximal particle number. The value 'Zero' signifies no 
constraint, so that obviously at most one of these can larger than zero, both zero standing 
for no constraint at all. Other options control the amount of output produced. 

There are default input values given specifying QED(1+1) with massive fermions and 
a strong coupling (g = 2.38m f). If the inputs are internally consistent they are exchanged 
against the defaults (check this, if there are problems). A convenient deviation made in 
the programming code is that the momenta and correspondingly the harmonic resolution 
are stored and assigned twice the value as in the formulae to make them integers instead 
of half- integers. For example, a harmonic resolution 10 is in fact K = 20 in the code and 
input files. 

A standard test run should take about 0.2s CPU time. Matrices of dimension of 
about 1500 take no more than 2000s, but use about 10 Mbyte RAM. This is because the 
whole matrix is handled as an array, the simplest possibility for the small matrices we 
are dealing with. The CPU time needed to generate the Fock space is generally negligible 
(Test Run: 0.008s) in comparison to the time needed to calculate the LC Hamiltonian 
matrix elements (0.098s) and diagonalize the matrix (0.062s). The diagonalization time is 
cubic in the matrix dimension and clearly dominates in the consumption of CPU time. 

3.2 Program structure 

The main program is 'ms_main.f. It 

• sets the maximum value for the matrix dimension (= number of basis Fock states), 
number of partons, resolution and the number of lines on one output page by in- 
cluding 'ms_parameters.f. Default values are dimension 700, particle number 10, 
resolution 600 and pagelength 60 lines. 

• sets input data defaults. The values correspond to QED(1+1) with Periodic Bound- 
ary Conditions (BC = true), fermion mass rrif = 1.0, a rather strong coupling 
g = 2.38 and a small resolution K = 20. Using the switch 'fixed particle number 
=0' and 'maximal particle number =0' we include the complete Fock space, e.g. all 
states with particle number up to the highest theoretically possible one. The output 
is limited to one page for each of the Fock state basis elements, mass eigenvalues 
and wave functions. 
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• reads the input data in subroutine 'input' (viz. screen output), tests it and exchanges 
it against the defaults if consistent. 

• calls the master subroutine 'master'. 

• writes out to the file 'masses. out', first a leading page with the input data (and the 
defaults if these were taken for the calculation), then the requested number of pages 
of Fock state basis, mass eigenvalues and wave functions with a header repeating the 
most important input data at the top of each page. 

• as a check writes out the actual input values to standard i/o. 
In the master subroutine 'master' 

• the Fock space basis is constructed using the subroutine 'con'. 

• the mass squared matrix parts 2KT, 2KC, 2KS and 2KF are calculated using the 
subroutines 'kinetic', 'contractions', 'seagulls' and 'forks'. 

• these parts are added to the correct mass squared matrix, and this is diagonalized 
using a NAGLIB routine. 

The important subroutines within 'con', 'kinetic', 'contractions', 'seagulls' and 'forks' are: 

• 'fermion' - constructs the charge singlet Fock states built from electrons and positrons. 
This essentially is the heart of the program and so we take the space here to ex- 
plain it in some detail. This subroutine is called from the routine 'con' and works 
as follows. For each fixed particle number allowed, the routine begins by building a 
reference state which consists of the lowest momentum electrons, namely a sequence 
of integers one after the other such that their sum does not exceed the total allowed 
momentum. This is stored as the first possible electron sub-state. From this state 
further states are now generated. The highest integer in the reference sequence is 
incremented. The total momentum of the new set is checked for total momentum 
saturation. If not saturated, the state is allowed and added to the electron sub-state 
array together with the sum of the particle momenta and the particle number. The 
procedure continues with incrementation of the highest electron momentum until 
total momentum is exhausted. The next-to-last momentum is then incremented, 
all higher momenta are reset to their (lowest) starting positions, and the routine 
proceeds with the highest momentum electron in an identical way to generate more 
states, always checking against total momentum saturation. Thus the program works 
its way down to the lowest momentum number. In this way an overly-large electron 
Fock space is built, but with all states satisfying the Pauli principle and momentum 
ordered. A positron space is built by copying the electron states. The two spaces are 
then joined by combining equal number electron and positron states. These com- 
bined states are then only checked against total momentum conservation as charge 
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conservation is already implicit. The stored particle numbers and total momenta of 
the sub-states allows for this to be very fast, so that further improvements were not 
implemented. 

• 'ep_graph' - calculates the matrix element contributions from electron-positron an- 
nihilation and electron-positron scattering graphs. It is called from the subroutine 
'seagulls'. It creates all possible matrix elements for these interaction operators by 
using the above-created space of states. In a second step, matrix elements are weeded 
out by checking first for particle number conservation. Then all possible participant 
combinations are tested for momentum conservation. Using the subroutine 'spec', 
the spectators are then checked for correct pairwise momentum conservation. As 
zero momentum transfers are unphysical in our theory, these are also excluded in 
a final step. The subroutines 'ee_graph' and 'pp_graph' are modelled on 'ep_graph' 
and all three taken together generate the 'seagull' matrix elements. 

• 'n_2_sec_e' - calculates the matrix element contributions from the particle number 
(n-2,n) sector graphs of an incoming electron undergoing pair production. It is called 
from the routine 'forks'. Again, as in 'ep_graph', all possible matrix elements are 
created and then non-allowed elements removed by checking that the number of 
'ingoing' and 'outgoing' particles differ by two and that momentum is conserved in 
the appropriate way. A further feature in this subroutine is the explicit inclusion 
of a relative minus sign between states with two indistinguishable electrons and/or 
positrons which arises directly from the Wick expansion. This subroutine is the 
template for the routines 'n_2_sec_p', 'njsec.eep' and 'n_sec_epp' and altogether give 
the 'fork' matrix elements of the light-cone Hamiltonian. 

• 'spec' - weeds out matrix elements for which spectators do not satisfy momentum 
conservation. The routine is given the participant configuration via the general par- 
ticipant vector 's' holding two ingoing and outcoming electron and positron momenta 
respectively. From that it is able to identify the remaining particles in a pair of states 
as spectators. For these pairwise momentum conservation is checked and when the 
check fails, the states are eliminated. 

• 'diag' - handles the calling of the general black box diagonalisation routine of the 
NAGLIB 'F02ABF'. We choose that routine because the QED(1+1) mass squared 
operator is real and symmetric (being hermitian) . It reduces the matrix to real sym- 
metric tridiagonal form by Householder's method and then uses the QL-Algorithm 
to extract eigenvalues and eigenvectors. A way of checking the results is to change 
this to routines 'F02AGF' and 'F02AXF', which deal with general nonsymmetric and 
hermitian matrixes, respectively. The matrix is stored as an array, which causes this 
method to use rather huge amounts of storage. However it is reliable, conceptionally 
simple and obtains the full spectra and wavefunctions. Note that the eigenvectors 
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are stored in columns in an array called 'vector'. 



3.3 Program package components 

In the package there are the following Fortran files: 

1. ms_main.f 

2. ms_parameters.f 

3. ms_master_sub.f 

4. ms_construct_sub.f 

5. ms_kinetic_sub.f 

6. ms_contractions_sub.f 

7. ms_seagulls_sub.f 

8. ms_forks_sub.f 

9. ms_spectator_sub.f 

10. ms_diagonalization_sub.f 

11. ms_in-out_sub.f 

12. ms_input_data 

File 1 contains the main program, file 2 the parameter definitions used by all routines and 
the files 3-11 the subroutines with self-explanatory file names. In file 12 the input of the 
test run is given in a form already prepared for automatized use. 



3.4 Reliability and Uses 

The routines shown here have already been applied to obtain mass spectra for various 
coupling constants and with different mass normalizations. The confirmation of their 
reliability comes from the very good comparison of the results with spectra from alternate 
analytic and numerical methods. In particular, these are []i"7| , 18, [l£|. Moreover we can 
directly compare against the m/ = case of the exact Schwinger model p0| . where the 
lowest mass eigenvalues should be 1 and 2. Even a resolution of 10 (K = 20) already 
yields values of mi = 0.99875 and 777-2 = 1.9994. The wavefunctions obtained with this 
method can be used to extract structure functions [15]. A novel example of uses for the 
so-obtained spectra is the computation of thermodynamic quantities for the given theory 
at finite temperature pT| . 

A generalization of this Abelian gauge theory model would be to include extra flavours. 
Moreover these routines are already the basis for application to other field theories in 1+1 
dimensions such as QED with massive photons |22]|. 



3.5 Test Run 

In Appendix (|A|) we give a sample screen output with the sequence of data inputting, 
checking and confirmation. The input values used here are the default values. The data 
file output of a test corresponding to this default input is shown in Appendix (Q). This 
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should be the basis of checking a test run of the program package. As a second option, in 
Appendix (|C]) we give a suggested input file which can be piped into the program while 
running our routines in batch mode. 

4 Concluding Remarks 

In this paper we have described a program to solve gauge theories for the bound state 
spectra and wavefunctions via the formal method of Discretized Light-Cone Quantiza- 
tion. The concrete test case for which we gave details for was quantum electrodynamics 
in one space and one time dimension. This model theory shares many features with 
the real problem to which DLCQ is eventually intended to be applied, namely quantum 
chromodynamics. These common features are the presence of Lorentz and gauge symme- 
try, fermions and a linear confining potential. However, the absence of spin, color and 
transverse momentum quantum numbers allow a simple exposition of the essentials of the 
programming procedure. We thus descibed routines which constructed a Fock space of 
electrons and positrons, and which then computed matrix elements of the Hamiltonian for 
charge singlet combinations of the Fock states. Diagonalization of this Hamiltonian has 
been used to generate, for example, spectra of bound states of electron-positron pairs in 
two-dimensions . 

Extension to non-Abelian gauge theories, namely QCD in two dimensions, is techni- 
cally more complicated in that color statistics must be introduced into the Fock space 
construction and more complicated interactions built into the code. However these repre- 
sent a 'grafting' onto the type of code we explicitly discussed here. Further extension to 
include more space-dimensions is also feasible in the same spirit pursued here, where color 
is compounded by the presence of genuinely dynamical gluons in the Fock space construc- 
tion. In both cases, the methods used here have been the basis for codes developed in, for 
example, pSfl . These two examples far from exhaust the variations of simplifications on 
full QCD which could be constructed in order to build toward the full problem and for 
which DLCQ, as a general method, could be insightful. As a further example we men- 
tion the novel application to two-dimensional Yang-Mills theory coupled to scalar adjoint 
matter by ||. For a treatment of full QCD, the major problems which could change some 
aspects of the programming method we describe are the renormalization and zero mode 
problems. These difficult formal problems are under active consideration. 
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l Screen output 

i 



The Output of 



on your terminal screen 



please give input data: 

Next input data set? (l=yes) 

boundary conditions 
(antiperiodic = T, periodic = F) 
boundary conditions are 



harmonic resolution 
harmonic resolution is 



20 



fermion mass (with decimal point) 
fermion mass is 



1.00 



coupling constant (with decimal point) 
coupling constant is 



2.38 



particle number limits 
(one of these must be zero; 
both zero is a full fock space request) 

fixed particle number 
fixed particle number is 

maximal particle number 
maximal particle number is 

front page (l=yes) 
front page printed 

maximal number of Fock space pages 
maximal number of Fock space pages is 

maximal number of eigenvalue pages 
maximal number of eigenvalue pages is 

maximal number of wavefunction pages 
maximal number of wavefunction pages is 



47 
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check output - you calculated with data: 



BC - boundary conditions: 



True (antiperiodic) 



res - harmonic resolution: 



20 



m_f - fermion mass: 



1.00 



c_c - coupling constant: 



2.38 



fixed particle number: 

maximal particle number: 6 

Output is: 1 of max 2 page(s) of Fock space 

1 of max 1 page(s) of mass values 

1 of max 12 page(s) of eigenvectors 

CPU-time for the DLCQ matrix. Dimension is 42 

total time: 0.168 s fock space: 0.008 

matrix: 0.098 s diagonalisation: 0.062 



The Output of 'ms' on your terminal screen 



please give input data: 



Next input data set? (l=yes) 



B Data output 



'ms' - output 



15-Jul-94 

2 BC: True 

res= 20 

3 

4 



m_f= 1.00 



c_c= 2.38 



18:07:06 
Nmax= 6 



+- 



'ms' 



program 'ms' calculated the mass spectrum of QED 1+1 
with the accepted input data: 



BC - boundary conditions: 



True 
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(antiperiodic) 
14 
I 

15 
I 

16 
I 

17 
I 

18 
I 

19 
I 

20 
I 

21 
I 

22 
I 

23 
I 

24 
I 

25 
I 

26 
I 

27 
I 

28 
I 

29 
I 

30 
I 

31 
I 

32 
I 

33 
I 

34 



res - harmonic resolution: 

m_f - fermion mass: 

c_c - coupling constant: 

fixed particle number: 

maximal particle number: 

full Fock space - max . part . number : 



20 



1.00 



2.38 



Output is: 1 of max 2 page(s) of Fock space 

1 of max 1 page(s) of mass values 
1 of max 12 page(s) of eigenvectors 

CPU-time for the DLCQ matrix. Dimension is 42 . 

total time: 0.168 s fock space: 0.008 s 

matrix: 0.098 s diagonalisation: 0.062 s 



page 2 : : : : 
'ms' - output 



61 

05-Jul-94 

62 BC: True 

res= 20 

63 

64 



m_f= 1.00 



c_c= 2.38 



18:07:04 
Nmax= 6 



65 
42 I 
66 



I Fock space basis used by 'ms' 



42 of 



67 



momenta 
69 



state mass~2 part lowest 5 electron momenta lowest 5 positron 



70 


1 


21.05 


2 




1 


19 












71 


2 


7.84 


2 




3 


17 












72 


3 


5.33 


2 




5 


15 












73 


4 


4.40 


2 




7 


13 












74 


5 


4.04 


2 




9 


11 












75 


6 


4.04 


2 




11 


9 












76 


7 


4.40 


2 




13 


7 












77 


8 


5.33 


2 




15 


5 












78 


9 


7.84 


2 




17 


3 












79 


10 


21.05 


2 




19 


1 












80 


11 


48.00 


4 


3 


1 


1 












81 


12 


34.87 


4 


3 


1 


3 












82 


13 


32.48 


4 


3 


1 


5 












83 


14 


31.75 


4 


3 


1 


7 












84 


15 


45.54 


4 


5 


1 


1 












85 


16 


32.48 


4 


5 


1 


3 












86 


17 


30.22 


4 


5 


1 


5 












87 


18 


44.68 


4 


7 


1 


1 












88 


19 


31.75 


4 


7 


1 


3 












89 


20 


29.71 


4 


7 


1 



15 
13 
11 

9 
13 
11 

9 
11 



14 



5 
90 

1 
91 

3 
92 

1 
93 

3 
94 

1 
95 

1 
96 

1 
97 

3 
98 

5 
99 

1 

100 
3 

101 
1 

102 
3 

103 
1 

104 
1 

105 
1 

106 
3 

107 
1 

108 
1 

109 
1 

110 
1 

111 
1 



21 


44 


44 


4 


9 


1 


22 


31 


75 


4 


9 


1 


23 


44 


68 


4 


11 


1 


24 


32 


48 


4 


11 


1 


25 


45 


54 


4 


13 


1 


26 


48 


00 


4 


15 


1 


27 


32 


48 


4 


5 


3 


28 


19 


56 


4 


5 


3 


29 


17 


52 


4 


5 


3 


30 


31 


75 


4 


7 


3 


31 


19 


05 


4 


7 


3 


32 


31 


75 


4 


9 


3 


33 


19 


56 


4 


9 


3 


34 


32 


48 


4 


11 


3 


35 


34 


87 


4 


13 


3 


36 


29 


71 


4 


7 


5 


37 


17 


52 


4 


7 


5 


38 


30 


22 


4 


9 


5 


39 


32 


48 


4 


11 


5 


40 


31 


75 


4 


9 


7 


41 


60 


19 


6 


5 3 


1 


42 


60 


19 


6 


7 3 


1 



: : page 3 : : : : 
'ms' - output 



121 
05-Jul-94 

122 BC: True 

res= 20 

123 

124 



m_f= 1.00 



c_c= 2.38 



18:07:04 
Nmax= 6 



7 
7 
5 
5 
3 
11 
9 
7 
9 
7 
7 
5 
5 
3 
7 
5 
5 
3 
3 
3 
3 



+ 

42 of 



125 
42 I 
126 



I Mass~2 eigenvalues calculated by 'ms' 



127 

128 
17.5337 

129 
22.8458 

130 
27.8514 

131 
38.3057 

132 
59.4818 

133 



181 
05-Jul-94 

182 BC: True 

res= 20 

183 

184 



1: 3.2217 7. 

9: 17.7096 17. 

17: 22.9126 23. 

25: 28.9435 29. 

33: 41.4285 43. 

41: 66.6438 75. 
: page 4 ::::::: 
ms' - output 



0677 
9037 
7122 
4743 
1637 
3417 



11 


2069 


11 


7232 


12 


3816 


14 


4305 


16 


0938 


19 


7839 


20 


0023 


21 


6207 


21 


7420 


21 


8782 


23 


9960 


25 


2102 


25 


7435 


26 


1450 


26 


6680 


33 


4945 


33 


6035 


35 


1072 


35 


3365 


37 


4794 


46 


1262 


46 


9805 


49 


7538 


51 


6405 


57 


6122 



m_f= 1.00 



c_c= 2.38 



18:07:04 
Nmax= 6 



185 

1 I 

186 



I 'ms' wavef unction (times 10) 



page 



1 of 



+ 

7 

16.09 



187 
188 
8 

189 
17.53 

190 

191 
1.6289 

192 
0.4102 

193 
-2.2542 

194 
0.6297 

195 
0.1992 

196 



state 
mass"2 



1 

3.22 



2 

7.07 



3 

11.21 



4 

11.72 



5 

12.38 



14.43 



1: 


1 


9568 


-2 


9886 


-3 


3115 





9786 


-1 


6094 


2 


6508 


-2 


7312 


2: 


2 


7551 


-3 


9533 


-3 


6808 





9977 


-1 


5350 


1 


8788 


-1 


4141 


3: 


3 


2898 


-3 


8657 


-1 


6452 





1565 





0947 


-2 


2260 


2 


7494 


4: 


3 


6375 


-2 


8099 


1 


6692 


-1 


1390 


2 


4319 


-1 


4712 


2 


6845 


5: 


3 


8105 


-1 


0277 


4 


1583 


-1 


7331 





4649 


-0 


6044 


-2 


4939 


6: 


3 


8105 


1 


0277 


4 


1583 


1 


7331 


-0 


4649 





6044 


-2 


4939 



15 



0.1992 








197 


7: 


3 


6375 


0.6297 








198 


8: 


3 


2898 


-2.2542 








199 


9: 


2 


7551 


0.4102 








200 


10: 


1 


9568 


1.6289 








201 


11: 





0508 


-0.7516 








202 


12: 





0107 


1.8602 








203 


13: 





0012 


-2.6779 








204 


14: 


-0 


0002 


3.8333 








205 


15: 





0484 


1 . 5882 








206 


16: 





0083 


-0.8839 








207 


17: 


-0 


0004 


. 9497 








208 


18: 





0275 


-0.5414 








209 


19: 


-0 


0003 


0.9056 








210 


20: 


-0 


0021 


-0.7469 








211 


21: 





0000 


0.0000 








212 


22: 


-0 


0092 


-0.0200 








213 


23: 


-0 


0275 


0.5414 








214 


24: 


-0 


0123 


-1.5719 








215 


25: 


-0 


0484 


-1.5882 








216 


26: 


-0 


0508 


0.7516 








217 


27: 





0123 


1.5719 








218 


28: 





0021 


-1.7284 








219 


29: 


-0 


0001 


2.1043 








220 


30: 





0092 


0.0200 








221 


31: 





0000 


0.0000 








222 


32: 





0003 


-0.9056 








223 


33: 


-0 


0021 


1 . 7284 








224 


34: 


-0 


0083 


0.8839 








225 


35: 


-0 


0107 


-1.8602 








226 


36: 





0021 


. 7469 








227 


37: 





0001 


-2.1043 








228 


38: 





0004 


-0.9497 








229 


39: 


-0 


0012 


2.6779 








230 


40: 





0002 


-3.8333 








231 


41: 





0000 


-0.0831 








232 


42: 





0000 


-0.0831 









233 



2 


8099 


1 


6692 


1 


1390 


3 


8657 


-1 


6452 


-0 


1565 


3 


9533 


-3 


6808 


-0 


9977 


2 


9886 


-3 


3115 


-0 


9786 





1227 





5743 


-0 


2382 





1923 





4471 


-0 


3046 





1634 





1320 


-0 


1184 





0651 


-0 


0075 


2 


4650 





3630 





9904 


-0 


5019 





3870 





4591 


-0 


4231 





2494 





0161 


2 


6933 





5813 





7688 


-0 


4837 





4728 





0473 


2 


4318 





1830 


-0 


1102 


-0 


1730 





6695 





0000 


2 


0339 





3973 


-0 


3850 


-0 


4348 





5813 


-0 


7688 


-0 


4837 





2058 


-0 


4299 


-0 


3354 





3630 


-0 


9904 


-0 


5019 





1227 


-0 


5743 


-0 


2382 





2058 





4299 


-0 


3354 





2121 





1534 


-0 


1939 





0895 


-0 


0056 


3 


8877 





3973 





3850 


-0 


4348 





3037 





0000 


3 


4922 





4728 


-0 


0473 


2 


4318 





2121 


-0 


1534 


-0 


1939 





3870 


-0 


4591 


-0 


4231 





1923 


-0 


4471 


-0 


3046 





1830 





1102 


-0 


1730 





0895 





0056 


3 


8877 





2494 


-0 


0161 


2 


6933 





1634 


-0 


1320 


-0 


1184 





0651 





0075 


2 


4650 


-0 


0001 





0045 





0106 





0001 





0045 


-0 


0106 



-2 


4319 


1 


4712 


2 


6845 


-0 


0947 


2 


2260 


2 


7494 


1 


5350 


-1 


8788 


-1 


4141 


1 


6094 


-2 


6508 


-2 


7312 





4703 


-1 


4469 


1 


1210 





7164 


2 


2283 


-0 


6816 


-2 


1230 


-1 


0377 


-0 


9995 


1 


2108 


-0 


1620 


1 


0982 


1 


0885 





9985 


-0 


1918 


-1 


6197 


1 


4899 


-2 


0400 


-1 


0611 


-1 


3940 





5794 


-1 


2287 





3499 


-2 


1919 


-0 


3335 


1 


4770 





1299 


-2 


6526 


-1 


3541 





8304 


-0 


1441 





3881 





0000 


-1 


7352 


1 


6020 


1 


7598 


-1 


2287 





3499 


2 


1919 





8182 


2 


7150 





6096 


1 


0885 





9985 





1918 





4703 


-1 


4469 


-1 


1210 





8182 


2 


7150 


-0 


6096 


-3 


1555 


-1 


6260 


-1 


3350 


1 


9939 


-0 


2040 





5719 


-1 


7352 


1 


6020 


-1 


7598 


-1 


4084 


-1 


8976 





0000 


-0 


3335 


1 


4770 


-0 


1299 


-3 


1555 


-1 


6260 


1 


3350 


-1 


6197 


1 


4899 


2 


0400 





7164 


2 


2283 





6816 


-2 


6526 


-1 


3541 


-0 


8304 


1 


9939 


-0 


2040 


-0 


5719 


-1 


0611 


-1 


3940 


-0 


5794 


-2 


1230 


-1 


0377 





9995 


1 


2108 


-0 


1620 


-1 


0982 





0180 





0109 


-0 


0823 


-0 


0180 


-0 


0109 


-0 


0823 



C Batch job input data 



1 1 ! input file for 'ms' 

2 .true. !BC boundary condition 

3 20 Ires harmonic resolution 

4 1.0 !m_f fermion mass 

5 2.38 !g coupling constant 

6 !N_fix fixed parton number 

7 !N_max maximal parton number 

8 1 !pages(l) front page print option 

9 1 !pages(2) Fock space print option 

10 1 Ipages (3) eigenvalue print option 

11 1 !pages(4) wavefunction print option 

12 -1 ! input file for 'ms' 
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